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Atmosphere, Ocean, and Land: 
Critical Qaps in 'Earth System Models 

Ronald G. Prinn and Dana Hartley 


Introduction 

The earth’s unique environment for life is determined by an interac- 
tive system comprising the atmosphere, ocean, land, and the living 
organisms themselves. We know from observations that this system is 
not static but changing. Our ability to understand past changes and 
predict future changes is hampered by insufficient understanding of at 
least five fundamental processes in the earth system: (1) convection, 
condensation nuclei, and cloud formation; (2) oceanic circulation and 
its coupling to the atmosphere and cryosphere; (3) land surface hydrol- 
ogy and hydrology-vegetation coupling; (4) biogeochemistry of green- 
house gases; and (5) upper atmospheric chemistry and circulation. 

We need to improve our knowledge of each of these processes 
through a combination of observations and theory (model develop- 
ment). We then need to incorporate realistic models for these 
processes into climate prediction models with special attention to 
potential feedbacks to climate change involving these processes (Fig- 
ure 1). In each case It is not just a matter of explaining the present 
state but also of understanding how each of the processes con- 
tributing to this state will change when one or more of the external 
forcing mechanisms changes. In this paper we briefly review current 
knowledge and pinpoint some of the major areas of uncertainty for 
each of the above five fundamental processes. 


Convection, Condensation Nuclei, and Cloud Formation 


Clouds and water vapor play a major role in the heat budget of 
the earth’s surface. Clouds contribute to surface cooling by reflect- 
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ing sunlight back to space. In fact, cloud albedo doubles the albedo 
of the entire earth over what it would be in the absence of clouds 
(Ramanathan et al., 1989). Furthermore, both water vapor and 
clouds contribute to surface heating by absorbing and re-emitting 
infrared radiation (greenhouse effect). Indeed, the contribution to 
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Figure 1 . Schematic illustrating required research Jor improving current climate models. 
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the greenhouse effect of clouds and water vapor is far greater than 
that caused by C0 2 (see Table 1 for chemical names and formuli) 
and the other long-lived greenhouse gases. 


Table 1 : Chemicals named in this chapter 


12C 

carbon - 12 

13 C 

carbon- 13 

CFCL 3 

trichlorofiuoromethane (CFC-11) 

CF 2 C 12 

dichloro-difluoromethane (CFC-12) 

CF 2 C1CFC1 2 

l-chlorodifluoro-2-dichlorofluoroethane (CFC-1 13) 

chcif 2 

chlorodifluoromethane (HCFC-22) 

chci 2 cf 3 

1 -dichloro-2-trifluoroe thane (HCFC- 123) 

ch 2 fcf 3 

1 -trifluoro-2-fluoroethane (HCFC- 134a) 

ch 3 cci 3 

I-trichloroethane (methyl chloroform) 

CH 3 5CI 3 (DMS) 

dimethyl sulfide 

CH 4 

methane 

Cl 

chlorine atom 

CIO 

chlorine monoxide 

ciono 2 

chlorine nitrate 

ClOOCl 

chlorine monoxide dimer 

Cl 2 

molecular chlorine 

CO 

carbon monoxide 

co 2 

carbon dioxide 

HC1 

hydrogen chloride 

HN0 3 

nitric acid 

ho 2 

hydroperoxyl free radical 

h 2 o 

water 

h 2 o 2 

hydrogen peroxide 

NO 

nitric oxide 

no 2 

nitrogen dioxide 

n 2 o 

nitrous oxide 

0 { 1 D) 

electronically excited oxygen atom 

OH 

hydroxyl free radical 

°2 

molecular oxygen 

°3 

ozone 

RH 

generic hydrocarbon 

ro 2 

generic organic peroxy free radical 

ROOH 

generic organic hydroperoxide 

so 2 

sulfur dioxide 


The net effect of clouds is determined by the opposing effect of the 
albedo of the cloud and its opacity to infrared radiation. Theory and 
observations have shown that the net effect of clouds is to lower the 
surface temperature relative to the cloud-free state (Somerville and 
Remer, 1984; Ramanathan et al., 1989). However, there are still dif- 
ficulties in representing the present cloud forcing in large-scale gen- 
eral circulation models (GCMs). A recent comparison between GCM 
results and present-day observations illustrates some of the weak- 
nesses in present parameterizations of cloud formation (Kiehl and 
Ramanathan, 1990). Figure 2 in particular shows that there is too 
much solar radiance reaching the surface near the poles and too 
little near the tropics in the community climate model (CCM) of the 
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Figure 2. Shortwave cloud radiative forcing (SWCF) as a Junction of latitude in 
the western North Atlantic from ERBE data and from the NCAR community 
climate model (adapted from Kiehl and Ramanathan, 1990). 


National Center for Atmospheric Research as compared to Earth 
Radiation Budget Experiment (ERBE) data. Some of this difference 
is associated with the surface albedo’s being too large in the model 
tropics and some with obvious gaps in the ability of the CCM to 
parameterize the cloud formation process. Furthermore, for predic 
tions of future climate change, we want to be capable of determining 
how the cloud /radiative balance may change, which requires a fun- 
damental understanding of the cloud formation process. We then 
need realistic models for linking cloud microphysics and atmos- 
pheric dynamics to yield accurate predictions of the global three- 
dimensional structure of clouds. 

The microphysics of clouds is determined to a large extent by the 
ambient concentrations of cloud condensation nuclei (CCN). Over 
land the presence of CCN is often due to dust or aerosol emissions by 
humans. Over the oceans the predominant source of CCN is either 
land-derived CCN and CCN precursors like S0 2 advected off the 
land, or the CCN precursor dimethyl sulfide (DMS), which is pro- 
duced in the surface ocean itself by planktonic algae (see, e.g. , Charl- 
son et al., 1987). The more CCN available, the smaller the cloud 
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droplets will be. Twomey (1977) demonstrates in addition that for a 
fixed volume of cloud water, the smaller the particles, the greater the 
albedo. Charlson et al. (1987) have specifically calculated (Figure 3) 
the change in cloud-top albedo as compared to a reference cloud as 
the droplets either grow or become smaller with the same total liquid 
water content. The top axis in Figure 3 indicates how this change is 
related to the number of droplets. Thus, as the number of droplets 


Relative number density of cloud droplets N/Nq 



Figure 3. Contours showing change of visible albedo caused by changing 
droplet number density N or equivalently droplet effective radius r e jj 
relative to a reference cloud (Nq, r Ijj). Albedo at top of reference cloud 
given on left-hand vertical scale (adapted from Charlson et al, 1987 ; 
reprinted by permission from Nature, © 1987 Macmillan Magazines Ltd.). 
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increases, as would be the case if more CCN were available due to 
increased DMS emissions, the albedo increases. Therefore, any cli- 
mate changes that influence the emissions of DMS over the oceans 
may induce a cloud-climate feedback associated with them. Under- 
standing how CCN may change is therefore an important aspect of 
predicting future cloud cover and rainfall. 

The three-dimensional structure of optically thick clouds is 
important. For very deep clouds, the downward emission from the 
warm cloud base significantly exceeds the upward emission from 
the cool cloud top. Very shallow clouds, on the other hand, exhibit 
much smaller differences between these fluxes, so the greenhouse 
effect due to them is much smaller. This influences the competition 
between the albedo and greenhouse effect of clouds. For example, 
Ramanathan and Collins (1991) have recently argued that relatively 
thin, highly reflective cirrus clouds formed during El Nino may effec- 
tively limit the magnitude of the sea surface temperature increase 
being forced by the greenhouse effect of water vapor and clouds. The 
three-dimensional structure of clouds depends on both mesoscale 
processes (convection) and the larger-scale circulation and is poorly 
predicted in current climate models. 

Convection also has a large influence on the surface heat budget. 
Evaporation and subsequent convection serve to cool and dry the 
surface and lower troposphere. In turn, the middle and upper tro- 
posphere are warmed by condensation. The distribution of water 
vapor and also many chemically important trace species is strongly 
influenced by convection and its subsequent motions. As shown by 
Dickerson et al. (1987), the convective motions tend to transport 
tracers from the lower troposphere to upper levels. Figure 4 illus- 
trates the redistribution of carbon monoxide. To the left of the cloud 
the normal vertical distribution is evident, with the CO mixing ratio 
decreasing with height. However, within the cloud, the situation has 
reversed as the high-CO air in the boundary layer has been trans- 
ported to the upper troposphere, leaving the lower levels depleted in 
CO relative to the normal situation. 

For climate prediction models we clearly need realistic models for 
convection, Including its roles in heat, momentum, and trace 
species (e.g., water vapor) transport, and in cloud formation. The 
first step in developing a cloud model is to understand the basic 
processes, and the second step is to parameterize these processes 
for inclusion in GCMs, in which they are subgrid-scale phenomena. 
Previous attempts at producing convective parameterizations have 
had their drawbacks. For example, Manabe et al. (1965) used a 
straightforward convective adjustment scheme that adjusts the 
atmospheric lapse rate back to an observed neutral state for satu- 
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Figure 4. Concentration of carbon monoxide (ppbv) as a function of altitude 
and latitude over Oklahoma City, Oklahoma, both outside and within a 
thunderstorm (adapted from Dickerson et al„ 1987; © 1987 by theAAAS). 


rated regions. Unfortunately, this means instability can build up in 
unsaturated regions. A recent parameterization by Betts (1986) 
allows for adjustment in unsaturated regions, but it still does not 
predict the profiles of temperature or water vapor. Instead, these are 
adjusted to observed profiles. Emanuel (1990) allows for the predic- 
tion of water vapor content but relies on a specified vertical distribu- 
tion of precipitation. The major drawback of the above schemes is 
that they are all based on adjusting the atmosphere back to an 
observed structure. For purposes of climate change models, this 
type of parameterization is quite suspect since certain potentially 
important climate changes are omitted. The Kuo (1974) scheme is 
more physical. In this scheme convection occurs when there is 
moisture convergence. However, this also has its problems. Regions 
can get extremely unstable, yet unless there Is a transport of mois- 
ture into that region, no convection occurs. Realistic models for con- 
vection and cloud formation are clearly an important requirement 
for GCMs to be used for future climate change predictions. 
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Oceanic Circulation and Its Coupling to the 
Atmosphere and Cryosphere 

The oceanic circulation has major influences on the climate system. 
In the energy balance of the earth, the ocean circulation is important 
in the heat transport from equator to pole. Figure 5 shows the heat 
• transport for both the oceans and the atmosphere as determined by 

s Carissimo et al. (1985) and earlier workers. The oceanic and atmos- 

pheric heat transport have comparable strengths. Thus, it is evident 
| that the ocean is crucial to consider and without it the latitudinal tem- 

l perature gradient would probably be significantly larger. The ocean cir- 

culation also contributes to the determination in space and time of the 
I latent, sensible, and radiative heating of the troposphere. For studying 

j possible climate changes, the changes in water vapor input from the 

I ocean with changing sea surface temperature are a critical feedback 

s (Ramanathan, 1981). Another way the ocean circulation affects the 

heating in the troposphere is in its interaction with sea ice. Figure 6 
gives some results from Hibler (1990) and earlier papers which show 
how important it is to consider the ice dynamics. There is about 50% 
more ice when there is no ice dynamics. The greater the extent of the 
sea ice, the more it can decrease the heat flux between the ocean and 
the atmosphere (Bentley, 1984). Since generally, in the polar regions, 
the ocean is warmer than the atmosphere, the extent of sea ice can 
greatly influence the temperature of the tropospheric air. Furthermore, 
the albedo of sea ice is about 0.9 whereas that of sea water is about 
0.1. Thus, the radiative heating of the troposphere is influenced by sea 
ice cover, which in turn is affected by the ocean circulation. 

The ocean circulation is essential for transporting nutrients to the surface 
of the open oceans in order to feed phytoplankton and other marine oigan- 
isms (Broecker and Peng, 1982). An example of the importance of this process 
is seen during an El Nino year when the usual East Pacific equatorial 
upwelling is suppressed. As a consequence there are no new nutrients, and 
the fish die off, which is devastating for the Peruvian fishing industry (Philan- 
der and Rasmusson, 1985). 

Finally, the ocean serves as a sink for thermal energy and carbon 
dioxide as they increase in the atmosphere. The ability of the ocean 
to do this depends greatly on the strength of its circulation. Figure 7 
is a result from Hansen et al. (1985) showing the different tempera- 
ture increases expected with different vertical diffusion rates in a 
highly simplified ocean model. Obviously, the ocean serves as a 
large heat sink since its heat capacity is so much greater than the 
atmosphere’s, but its ability to do so depends on its circulation. This 
is also true in the uptake of C0 2 . Siegenthaler (1983) explored the 
effects of circulation on the uptake of C0 2 , and one of his results is 
shown in Figure 8 which illustrates how much more C0 2 is taken 


‘RpnaCd g. Trinn and Dana Manky 


17 



Latitude 


Figure 5. Annual meridional transport of energy (10 15 W) in the atmos- 
phere plus ocean (T A + T ( f, atmosphere (T A ). and ocean (T 0 ) computed 
by Cartssimo et al. (1985) and compared to previous investigators: (+) 
Oort and Yonder Hoar, 1976; (x) Trenbenh, 1979; (O) Newell et al, 1972 
(adapted from Carissimo et aL, 1985; ©American Meteorological Society). 




Average ‘Basin Ice ‘Thickness (m) 


18 


Modeling the ‘Earth System 



Figure 6. Comparison of Arctic basin average ice thickness (m)for a two- 
level model without dynamics and with dynamics (standard) (adapted 
from Hibler, 1 990; © American Meteorological Society). 


into the ocean when there is efficient downward advection in "out- 
crop” regions like the Norwegian and Weddell seas. 

The atmosphere also influences the ocean since the upper ocean 
circulation is forced by atmospheric winds. This coupling is exempli- 
fied by the El Nino/Southem Oscillation (ENSO) phenomenon, dur- 
ing which the warming of the eastern equatorial Pacific is accompa- 
nied by slowing of the trade winds and by an eastward extension of 
the western Pacific convection zone. Slowing of the trade winds 
decreases the upwelling associated with Ekman drift, thus con- 
tributing to the warming of the eastern Pacific. The assignment of 
cause and effect in this coupled ocean -atmosphere phenomenon is 
difficult to make (Philander and Rasmusson, 1985). 

The ocean circulation and sea ice are also coupled. When ice is 
formed, the ocean becomes more saline, tending to increase the 
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oceanic overturning (Hibler, 1990). If the ice is transported by the 
combined effects of the winds and the ocean circulation (which are 
related) and then melts in another region, this latter region becomes 
less saline, and overturning is inhibited there. One possible conse- 
quence of this could be changes in the northern Atlantic climate. 
The Gulf Stream transports warmer water northward. When it 
reaches the North Atlantic the water cools, thereby warming the tro- 
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Figure 7 . Ocean surface response time (time to reach 1 - e 1 of equilibrium 
response) for a one-dimensional diffusive ocean model as a function of 
climate sensitivity to doubled C0 2 ^ 2 * 00 ^/ or equivalently as a 

function of the climate feedback factor f Mixed layer depth is 100 m and 
results for various vertical diffusion coefficients k (cm 2 /s) are shown 
(adapted from Hansen et at, 1985; © 1985 by the AAAS). 
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Figure 8. Fractions of the C0 2 added exponentially with an e folding time of 22.5 
years which are taken up by the ocean (diffusion into mixed layer plus simu- 
lated advection into deep water via an outcrop) as a function of the relative area 
of the outcrop. Circles denote fractions of bomb-produced 14 C taken up by the 
different media for relative outcrop areas of 0.0 and 0. 1 (adaptedfrom Siegen- 
thaler, 1983). 

posphere. As the water cools, it sinks, to be replaced by more warm 
water from the south (Broecker, 1987; Broecker and Denton, 1990). 
This northward flow of warm water in a “conveyor belt” keeps 
Europe warm during its winters. However, if sea ice invaded the 
North Atlantic and melted there, this would inhibit sinking (e.g., by 
the same mechanism in which the thermohaline circulation weak- 
ened due to increased river runoff and precipitation in the coupled 
ocean-atmosphere climate model of Stouffer et al., 1989). This 
would tend to decelerate the conveyor belt and cool Europe. This is 
just a hypothesis, but it could mean that one signal of greenhouse 
warming may actually be a cooling in some locations, which draws 
attention to the different results that could arise in a coupled ocean- 
atmosphere model but not in a model with prescribed sea surface 
temperatures or a horizontally static oceanic mixed layer. In fact, 
many ocean models that do have a circulation force the model to 
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observed climatology, which short-circuits many potential feedbacks 
when the climate is changing. We ultimately need realistic models 
that predict the changes in the ocean circulation as the climate 
changes and the resultant feedbacks to the climate and nutrient 
budgets (i.e., phytoplankton populations). The phytoplankton play a 
role in the uptake of 0O 2 as well as the release of DMS, which, as 
discussed earlier, may be an important influence on the formation 
and radiative properties of clouds. 

Land Surface Hydrology and 
Hydrology- Vegetation Coupling 

The relationship between vegetation, soils, water flow, and 
weather patterns is complex and fundamental to life on the conti- 
nents. The large-scale correlations between seasonal vegetation 
index, deep convective clouds (Warren et al., 1986), and surface 
albedo (Ramanathan et al., 1989) are quite evident. Vegetation is a 
major contributor to the absorption of sunlight over the land and, 
thus, to the local and global radiation budget. The albedo over vege- 
tation is less than half that over bare land (Chamey et al., 1975). 
Changing land vegetation cover can thus alter the net solar absorp- 
tion. In addition, vegetation pumps water from soils into the atmos- 
phere through evapotranspiration, thereby regulating the relative 
amounts of water in the atmosphere and in the soils. Thus, the veg- 
etation affects the local and downstream water and latent heat bud- 
get. Shukla and Mintz (1982) illustrate the control vegetation could 
have on the regional climate. Figure 9 shows the significant effect 
soil moisture can have on precipitation patterns, with Figure 9a 
having a wet soil and Figure 9b having a dry soil boundary. Shukla 
and Mintz (1982) also show the effects of the soil moisture on sur- 
face temperature and pressure. The dry soil causes the land to be 
warmer and produces lower pressures because the energy from 
solar insolation cannot go into evaporation at the surface. Further- 
more, there are fewer clouds (since there is less water to feed them), 
and therefore even more solar radiation reaches the surface. Thus 
vegetation is a major influence on regional climate. 

On the other hand, vegetation is itself dependent on the climate. 
Different vegetation types thrive in different climates (Running and 
Coughlin, 1988; Clark, 1990, 1991). Some species are very sensitive 
to precipitation and solar input. Therefore, as these change, so will 
the vegetation. Furthermore, as regions become more diy the fire 
potential increases, and large areas of vegetation can potentially be 
cleared by burning. This can then influence what new vegetation will 
grow (Clark, 1990, 1991). Thus, changes in climate can both directly 
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Figure 9. Model simulated July precipitation (mm/ day) in (a) the 
wet soil and (b) the dry soil case (adapted from Shukla and 
Mintz, 1982 ; © 1982 by theAAAS). 

and indirectly alter the vegetation. One feedback mechanism is illus- 
trated in a well-known theory of desertification (Charney et al., 
1975). Charney et al. (1975, 1977) show that as vegetation is cleared, 
the albedo increases, evapotranspiration decreases, subsidence 
occurs, and the normal precipitation is suppressed. Thus, this feed- 
back dries out the surface, and recovery of the vegetation is severely 
Inhibited. It is apparent, therefore, that accurate models for predict- 
ing future climate change cannot merely specify the vegetation at the 
surface, but instead must include vegetation as an internal variable. 

For the current global GCMs, the effects of vegetation need to be 
parameterized since the scale of the relevant ecosystems is much 
smaller than that of a model grid box. These parameterizations must 
take account of local topography, water sources, and soil types since 
these are some of the crucial influences on vegetation cover. Various 
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simple vegetation /hydrology schemes in the GCMs of today give pre- 
dictions that do not agree well with observations or with one 
another on a regional scale. Figure 10 shows the different predic- 
tions of soil moisture from several different GCMs. For these calcu- 
lations, both the subgrid -scale convective and vegetation parameter- 
izations contribute to the distinct differences. 



Figure JO. Present-day simulations of soil moisture in winter for 
North America obtained in general circulation climate models 
from the National Center for Atmospheric Research , Geophysical 
Fluid Dynamics Laboratory , Goddard Institute for Space Stud- 
ies, U.K. Meteorological Office , and Oregon State University. 
Units are cm of water (adapted from Kellogg and Zhao, 1 988). 
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Biogeochemistry of the Long-Lived Greenhouse Gases 

Current concerns about future climate change are driven in large 
part by the undeniable observational evidence that several of the so- 
called long-lived greenhouse gases (C0 2 , CH 4 , N 2 0, chlorofluorocar- 
bons) are increasing presently at very significant rates. These 
increases, coupled with expected feedbacks increasing the concen- 
tration of the major greenhouse gas, the short-lived H 2 0 molecule, 
are predicted to lead to a significant warming over the next century 
(e.g., Dickinson and Cicerone, 1986). However, we lack the detailed 
biogeochemical and physical knowledge of individual sources and 
sinks for the long-lived greenhouse gases (with the exception of the 
chlorofluorocarbons) needed to explain quantitatively these present 
trends and to project them accurately into the future. 

The long-lived greenhouse gases can be conveniently divided into 
three categories based on their chemical reactivity: 

• First category: Some of the long-lived species are very inert (life- 
times exceeding 70 years) because they can be destroyed only by 
mixing up into the stratosphere above 25 kilometers altitude. 
Here, they can be dissociated by solar ultraviolet photons with 
wavelengths less than 250 nm. In this category are the industrial 
chlorofluorocarbons CFC1 3 , CF 2 C1 2 , and CF 2 C1CFC1 2 , which are 
widely used as refrigerants, propellants, plastic foaming agents, 
and solvents and which have lifetimes of around 45, 125, and 1 10 
years respectively at the present time (Golombeck and Prinn, 
1986, 1989, 1991, personal communication). Measurements 
made in the Atmospheric Lifetime Experiment/Global Atmos- 
pheric Gases Experiment (ALE /GAGE) network of automated sta- 
tions (located in Ireland, Oregon, Barbados, Samoa, and Tasma- 
nia; Prinn, 1988; Prinn et al., 1983) indicate that CFC1 3 , CF 2 C1 2 , 
and CF 2 C1CFC1 2 are currently increasing at rates of about 5%, 
5%, and 10% per year respectively over the globe (Figure 1 1). Also 
in this class is N 2 0, which has a present-day lifetime of about 
130 years and is increasing at 0.2 to 0.3% per year (Figure 11) 
and has both natural (soils, ocean) and anthropogenic (disturbed 
and fertilized soil, combustion) sources. Analysis of ALE/GAGE 
N 2 0 data using inverse methods yields the conclusion that tropi- 
cal sources of this gas compose 52-68% of the global total (Prinn 
et al., 1990). We need to understand these sources better and dis- 
cern their sensitivity to climate changes. Another point worth not- 
ing is that all of the gases in this first category, when dissociated 
in the stratosphere, yield reactive species that catalytically 
destroy ozone, and some of them play a major role in the 
processes causing the Antarctic ozone hole. Destruction of 
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Figure 11. Monthly average concentrations of various long-lived trace gases 
measured at the Cape Grim, Tasmania, station of ALE/GAGE. Left vertical 
axis gives numerical values for mixing ratios of all compounds except CH 4 , for 
which values are on the right vertical axis. Units are given on the data lines 
(from R. Prinn, P. Fraser, P. Stmmonds, R. Rasmussen, F. Alyea, D. Cunnold, 
A. Crawford, and R. Rosen, personal communication). 

stratospheric ozone in turn Increases ultraviolet dosages at the 
surface, affecting living organisms and influencing the oxidation 
process in the lower atmosphere by accelerating photodissocia- 
tion of 0 3 to produce O^D). 

• Second category: Long-lived species in this category are less inert 
(lifetimes typically less than 15 years) because they can be 
destroyed by reaction with OH in the lower atmosphere. In this 
category is atmospheric CH 4 with a lifetime of 9.6 years (Prinn et 
al., 1987), which is increasing today at about 1% per year (Figure 
11; see also Blake and Rowland, 1988). Methane (like N 2 0) has 
both natural (wetlands, tundra) and human-controlled (cattle, 
combustion, rice agriculture) sources (Ehhalt, 1988; Cicerone and 
Oremland, 1988). Also in this category is the widely used indus- 
trial solvent and cleaning agent methyl chloroform (CH 3 CC1 3 ) with 
a lifetime of 6.3 years (Prinn et al., 1987), which is increasing 
globally at about 4% per year (Figure 1 1). Finally, in this category 
are the hydrochlorofluorocarbons such as CHC1F 2 (16-year pres- 
ent-day lifetime), currently increasing at about 10% per year, and 
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CHC1 2 CF 3 (1.5 year steady-state lifetime), and CH 2 FCF 3 (15-year 
steady-state lifetime), proposed as industrial replacements for the 
chlorofluorocarbons CFCI 3 and CF 2 C1 2 . 

As reviewed by Crutzen (1979), the concentrations of OH radicals 
are determined to a significant extent by short-lived species 
whose life cycles begin and end on local and regional scales (e.g. , 
NO, N0 2 , 0 3 , CO, hydrocarbons, etc., in Figure 12). Thus the 
global atmospheric chemistry of the long-lived species in the sec- 
ond category is linked in very important ways to the local and 
regional-scale atmospheric chemistry of the short-lived species 
that determine OH levels. One important implication of this link- 
age is that the positive trends in species like CH 4 may in part be 



ical reactions. Thin arrows denote pathways for chemical reactions, 
while thick arrows denote flaxes of species from sources or to sinks. 
In continental air, rich in nitrogen oxides, NO is a major reactant for 
recycling H0 2 to OH and photodissociation of the resultant N0 2 is a 
major source for ozone. The dashed line enclosing NO and N0 2 sig- 
nifies that both are produced from surface sources and lightning. 

In remote marine air, where nitrogen oxide levels are low, the 
major recycling of H0 2 to OH occurs through peroxide (H 2 0 2 , ROOH) 
formation and subsequent photodissociation. Reaction ofH0 2 with 
0 3 provides another OH recycling mechanism in both environments. 
The OH and N0 2 are removed by formation of HN0 3 , which is 
scavenged by clouds, indicated by the dotted area 
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due to a decrease in OH radicals and thus a decrease in CH 4 loss 
rates rather than exclusively to an increase in CH 4 emission 
rates. Another implication is that the destruction of long-lived 
gases in this second category is weighted heavily to the hot tropi- 
cal lower atmosphere whose local and regional chemistry we are 
only beginning to be informed about (Donahue and Prinn, 1990). 

• Third category: C0 2 requires a category of its own. It is essentially 
chemically inert in the atmosphere and is cycled between atmos- 
phere, upper ocean, green vegetation, and soils on decadal time 
scales and between the upper ocean and deep ocean on a time 
scale of a few hundred years. It is currently Increasing at a rate of 
about 0.34% per year (Figure 13; Keeling, personal communica- 
tion) and possesses a significant annual cycle driven by uptake of 
C0 2 by photosynthesizing vegetation in spring and summer and 
release of C0 2 by decaying vegetation in fall and winter. Classical 
ideas, as reviewed by Oeschger and Siegenthaler (1988), had the 
oceans playing the major role as a sink for carbon dioxide. A typi- 
cal global budget has the current fossil fuel combustion source 
(5300 Tg C/yr) opposed by an oceanic sink (2300 Tg C/yr) to yield 
the observed atmospheric increase (3000 Tg C/yr). The land bios- 
phere was generally considered to be a small contributor to the 
budget, although estimates of the deforestation source ranged 
from 400 to 2600 Tg C/yr. A recent comparison of the observed 
C0 2 latitudinal gradient with a global model containing only an 
ocean sink (Tans et al., 1990) showed poor agreement (Figure 14). 
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Figure 13. Concentration of C0 2 (ppmv) measured at Mauna Loa 
Observatory, Hawaii (C.D. Keeling, personal communication). 
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Sine of Latitude 

Figure 14. Observed C0 2 mixing ratios (circles with 1 sigma error bars), a 
cubic polynomial fit to observed C0 2 (line a), and predicted C0 2 in three 
GCM model runs (lines b, c, d). In the model runs the ocean is the only C0 2 
sink, and this sink is computed in different ways in the three runs. A good 
ft (not shown) to the observations is obtained if a large Northern Hemi- 
sphere continental sink is included and the ocean sink decreased appropri- 
ately (adapted from Tans et al., 1990; © 1990 by the AAAS). 

A new, significant Northern Hemispheric sink is needed to fit the 
data, but the available ocean uptake estimates do not allow the 
Northern Hemispheric oceans to provide this sink. Tans et al. con- 
cluded therefore that a significant Northern Hemispheric continen- 
tal sink (about 2700 Tg C/yr) augmented by a smaller oceanic sink 
(about 1000 Tg C/yr) was needed to balance the atmospheric 
accumulation (3000 Tg C/yr), the fossil fuel input (5300 Tg C/yr), 
and the deforestation input (1400 Tg C/yr). Thus the ocean sink 
and the net biospheric sink (2700 -1400 = 1300 Tg C/yr) are com- 
parable. If correct, this implies that continental forest regrowth 
and soil carbon uptake will need to be understood and appropri- 
ately modeled. Of particular interest is the possibility of CO 2 fertil- 
ization. Perhaps related to this is the increasing amplitude of the 
C0 2 annual cycle in recent years (Figure 13). It is worth noting, 
however, that the available i^C/^C data do not indicate a bios- 
pheric sink of the size needed here (Keeling et al., 1989). 
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Better knowledge of the controlling mechanisms and magnitudes 
of the individual sources and sinks for the greenhouse gases is 
essential for developing the needed models for reliable prediction of 
future emissions and for inclusion of potential feedbacks between 
these source and sink rates and climate change. Of particular 
importance is improved understanding of the role of the biota as 
trace gas sources and sinks and the effects of climate change on 
these sources and sinks (see, e.g., Prinn, 1992). The potential for 
such biogeochemical feedbacks is seen in the analyses of gases 
trapped in ice cores, which show this effect over the past 150,000 
years where variations in concentrations of both methane and car- 
bon dioxide are correlated with variations of temperature (Chappel- 
laz et al., 1990). In particular, CH 4 and CO 2 are generally low during 
the ice ages and high during the warm Interglacial epochs (Figure 
15). This suggests the existence of a biogeochemical-climate feed- 
back whereby increases in temperature or greenhouse gas concen- 
trations lead to increases in greenhouse gas emissions or global 
warming respectively. 


Upper Atmospheric Chemistry and Circulation 

Changes in surface trace gas emissions, stratospheric-tropos- 
pheric exchange rates, and stratospheric temperatures and circula- 
tion will all affect the ozone layer. N 2 0, methane, industrial chloro- 
fluorocarbons like CFCI 3 and CF 2 Cl 2 , and water vapor are 
transported from the surface up into the stratosphere where photo- 
chemical reactions driven by solar ultraviolet radiation can decom- 
pose them to form a number of very reactive species (NO, N0 2 , Cl, 
CIO, OH, H0 2 ) which can catalytically destroy ozone. Because 
stratospheric ozone is produced essentially exclusively by ultraviolet 
dissociation of 0 2 , its rate of global production is set externally by 
the supply of ultraviolet radiation from the sun. The global rate of 
destruction of ozone depends, however, on the supply of catalysts, 
which is rapidly rising due to continued chlorofluorocarbon emis- 
sions at the surface. Removal of the catalysts from the stratosphere 
occurs through their (partially reversible) conversion to harmless 
reservoir species (HNO 3 , C10N0 2 , HC1, etc.), many of which can be 
transported down again to the troposphere ultimately to rain out or 
deposit at the surface. 

A dramatic illustration of the fact that the global chemical system 
is changing is provided by the Antarctic ozone hole phenomenon 
(Figure 16). This phenomenon can occur because temperatures in 



2000 



Figure 15 . Vostok ice-core analyses of CH 4 (upper curves, 2 sigma errors), temper 
ature relative to present as deduced from H 2 0 isotopes ( middle curve), and C0 2 
(louder curve , 2 sigma errors). Depth scale applies to C0 2 cm d CH 4 . Due to down- 
ward diffusion the gases are younger than the ice in which they are trapped 
(adapted from Chappellaz et aL, 1990 ; © 1990, Macmillan Magazines Ltd.). 
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Figure 16. October monthly mean total ozone measurements over Halley Bay 
as measured by Halley Bay Dobson and Nimbus 4 BUV (backscatter 
ultraviolet) instruments. The polar minimum values for TOMS (total ozone 
mapping spectrometer) for October are also shown. Dobson unit is 2.7 x 10 16 
ozone molecule per cm 2 (Schoeberl et al, 1989; Stolarski , 1988). 
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ANTARCTIC OZONE HOLE 



Figure 1 7. Schematic showing the vortex circulation, ice clouds, chlorine 
release reactions, and catalytic ozone destruction reactions involved in 
the Antarctic ozone hole. 


the meteorologically isolated air mass lying between 10 and 30 km 
altitude in the Antarctic winter become so low (less than 190 K) that 
ice crystals condense in this region (Figure 17). Decomposition of 
some of the catalyst reservoir species occurs on and within these ice 
crystals, producing CIO in the early Antarctic spring. The CIO levels 
become very high because HNO 3 remains in the ice and is therefore 
not available to contribute to the conversion of CIO back to its rela- 
tively harmless reservoir compound C 10 N 02 - The scarcity in the 
Antarctic spring of ultraviolet radiation capable of decomposing 0 3 
means that different catalytic cycles operate here than the ones 
operating in tropical regions, which require O atoms produced by 0 3 
photodissociation (e.g., Rowland and Molina, 1975). The reaction 
sequence shown in Figure 17 serves to efficiently destroy 0 3 without 
requiring 0 3 photodissociation and without destroying CIO (Molina 
etal., 1987). 

Theoretical models of the greenhouse effect indicate that a rise in 
the level of the greenhouse gases will tend not only to warm the 
earth’s surface and lower atmosphere but also to cool the stratos- 
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phere, due to either the rise in a greenhouse gas, as is the case for 
C0 2 , or the depletion of ozone, as in the case of rising chlorofluoro- 
carbons (e.g., Ramanathan et al., 1985; Wang et al„ 1991). This 
cooling of the stratosphere is expected to affect the ozone layer by 
decreasing ozone destruction in equatorial regions by slowing down 
the destruction reactions in these regions and increasing ozone 
destruction in the Antarctic (and also in the Arctic if cooling is suffi- 
ciently large) by increasing the extent of the polar ice clouds. These 
effects would be added to the increased ozone destruction expected 
due to the still rising chlorofluorocarbon concentrations. The 
stratospheric cooling by the greenhouse gases combined with these 
ozone layer changes will in turn affect the circulation at this level 
and conceivably at lower levels also. Also, changes in the intensity of 
tropical convection and midlatitude storms associated with a global 
warming are likely to alter significantly tropospheric-stratospheric 
exchange rates of trace gases involved in the ozone layer chemistry. 

We therefore need realistic models of stratospheric circulation 
and chemistry that can be coupled with models of the lower atmos- 
phere to improve the accuracy of tropospheric as well as stratos- 
pheric predictions. 


Concluding Remarks 

Understanding and modeling the above five fundamental 
processes provides a major interdisciplinary challenge with both 
observational and theoretical components. Several ongoing projects 
of the World Climate Research Program and the International 
Geosphere-Biosphere Program will all make large contributions. 
These projects include among others the International Global 
Atmospheric Chemistry (IGAC) project. Joint Global Ocean Flux 
Study (JGOFS), Global Energy and Water Cycle Experiment 
(GEWEX), Global Change and Terrestrial Ecosystems (GCTE) pro- 
ject, Tropical Oceans and Global Atmosphere (TOGA) project, and 
the World Ocean Circulation Experiment (WOCE). Later this decade 
the Earth Observing System (EOS) of the National Aeronautics and 
Space Administration will add a major satellite remote sensing capa- 
bility relevant to many of these projects. 
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